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We present a variational formulation for the Kardar-Parisi-Zhang (KPZ) equation that leads to a 
thermodynamic-like potential for the KPZ as well as for other related kinetic equations. For the 
KPZ case, with the knowledge of such a potential we prove some global shift invariance proper- 
0^ \ ties previously conjectured by other authors. We also show a few results about the form of the 

stationary probability distribution function for arbitrary dimensions. The procedure used for KPZ 
was extended in order to derive more general forms of such a functional leading to other nonlinear 
kinetic equations, as well as cases with density dependent surface tension. 
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£NJ , 1. Introduction 

' Phenomena far from equilibrium are ubiquitous in nature, including among many other, turbulence in fluids, interface 
and growth problems, chemical reactions, biological systems, as well as economical and sociological structures. During 
O ' the last decades the focus on statistical physics research has shifted towards the study of such systems. Among those 
studies, the understanding of growing kinetics at a microscopic as well as on a mesoscopic level constitutes a major 
challenge in physics and material science [Tong & Williams, 1994; Barabasi & Stanley, 1995; Halpin-Healy & Zhang, 
1995; Marsili et ai, 1996]. Some recent papers have shown how the methods and know-how from static critical 
phenomena have been exploited within nonequilibrium phenomena of growing interfaces, obtaining scaling properties, 
symmetries, morphology of pattern formation in a driven state, etc [Hentschel, 1994; Prahofer & Spohn, 2004; Lopez 
et ai, 2005; Fogedby, 2006; Ma et al, 2007; Castro et ai, 2007]. 

Even though it was (briefly) discussed in [Cross & Hohenberg, 1993], there is still a common belief that the nontrivial 
spatial-temporal behavior occurring in several nonequilibrium systems, originates from the non-potential (or non- 
variational) character of the dynamics, meaning that there is no Lyapunov functional for the dynamics. However, 
Graham and co-workers have shown in a series of papers [Graham, 1987; Graham & Tel, 1990] that a Lyapunov-like 
functional exists for a very general dynamical system, the complex Ginzburg-Landau equation. Such a functional 
\ is formally defined as the solution of a Hamilton-Jacobi-like equation, or obtained in a small gradient expansion 
- 1 ' [Descalzi & Graham, 1994]. The confusion associated with the qualification of nonvariational dynamics comes from 
the idea that the dynamics of systems having nontrivial attractors (limit cycle, chaotic) cannot be deduced from the 
. minimization of a potential playing the same role as the free energy in equilibrium systems [Cross & Hohenberg, 
t^J- ' 1993]. Nevertheless, this does not preclude the existence of a Lyapunov functional for the dynamics that will have 
local minima identifying the attractors of the system [Montagne et ai, 1996; Wio, 1997; Wio et ai, 2002; Wio & Deza, 
0\ ■ 2007]. However, once the system has reached an attractor that is not a fixed point, the dynamics proceeds inside the 
attractor driven by nonvariational contributions to the dynamical flow, that do not change the value of the Lyapunov 
■ functional, implying that the dynamics is not completely determined once the indicated functional is known. This 
situation has known examples even in equilibrium statistical mechanics [Hohenberg & Halperin, 1977]. Hence, the 
Lyapunov functional, or nonequilibrium potential (NEP) [Graham, 1987; Graham & Tel, 1990; Wio, 1997; Ao, 2004], 
plays the role in nonequilibrium situations of a thermodynamical-like potential characterizing the global properties 
of the dynamics: attractors, relative (or nonlinear) stability of these attractors, height of the barriers separating 
attractions basins, offering the possibility of studying transitions among the attractors due to the effect of (thermal) 
fluctuations. 

In a recent series of papers we have shown several results related to the obtention of the indicated NEP's for other 
system's classes: scalar and non-scalar reaction-diffusion systems [Bouzat & Wio, 1999; von Haeften et ai, 2004; 
von Haeften et at, 2005; von Haeften & Wio, 2007]. In particular we have exploited those results for the study of 
stochastic resonance [Gammaitoni et ai, 1998] in extended systems [Wio, 1997; Bouzat & Wio, 1999; von Haeften et 
ai, 2000; von Haeften et ai, 2004; von Haeften et ai, 2005; von Haeften & Wio, 2007; Tessone & Wio, 2007; Wio 
& Deza, 2007;]. In those works, we have analyzed problems of stochastic resonance in scalar and activator-inhibitor 
systems, in systems with local and nonlocal interactions, studied system-size stochastic resonance, etc. 
Here, and related to the kinetics of growing interfaces, we discuss the case of the Kardar-Parisi-Zhang equation (KPZ) 
[Kardar, Parisi & Zhang, 1986; Medina et ai, 1989]. This equation, that describes the evolution of h(x,t), a field that 
corresponds to the height of a fluctuating interface, reads 

= vV 2 h{x, t) + ± (Vh(x, t)) 2 + K + Z(x, t), (1) 
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where t) is a Gaussian white noise, of zero mean ((£(5, t)) = 0) and correlation t)£(x', t')) = 2e8(x — x')6(t — 
t'). As indicated above, this nonlinear differential equation describes fluctuations of a growing interface with a surface 
tension given by v, A is proportional to the average growth velocity and arises because the surface slope is parallel 
transported in such a growth process. 

Opposing to a claim in a recent paper [Fogedby, 2006]: The KPZ equation is in fact a genuine kinetic equation 
describing a nonequilibrium process in the sense that the drift v'S7 2 h + -|V/i • Vh — F cannot be derived from an 
effective free energy; in [Wio, 2007] it was shown that such a nonequilibrium thermodynamic-like potential 
(NETLP) for the KPZ equation exists. 

In this paper we present the approach to obtain the NETLP for the KPZ equation, and also show how it is possible to 
obtain such a NETLP for other related kinetic equations, as well as extend the procedure to other general situations. 
The organization is as follows. In the next Section we show how to obtain the NETLP for the KPZ case and, exploiting 
its knowledge, we discuss conjectures advanced in [Hentschel, 1994] and how they are fulfilled. Next we discuss about 
the form of the stationary (or asymptotic) probability distribution function, valid for any dimension. In the 
following Section we discuss the situation with a non-local interaction and symmetric kernels, showing how this form 
of extending the derivation procedure allows to obtain more general forms of kinetic equations. The following Section 
discusses the case of density dependent surface tension, while the last one presents some final comments. 

2. Kardar-Parisi-Zhang case 
2.1 Derivation 

In order to show how to obtain the above indicated functional, we start considering the following simple (general) 
scalar reaction-diffusion equation with multiplicative noise 

^<f>(x, t) = vV 2 (f>{x, t) + f(ct>(x, t)) + 4>{x, t)r){x, t), (2) 

where f(4>(x, t)) is a general nonlinear function. The noise r/(x, t) is a Gaussian white noise, of zero mean and intensity 
a. For this equation we assume the Stratonovich interpretation. 
It is known that the system in Eq. @ has the following NETLP 

x.t) ~\ 

f(v)dv+ V -{V^x,t)f\dx, (3) 
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where Vt indicates the integration range, and 

l^ <) = -^) + ^' i) ^ t); (4) 

where the contribution from the boundaries is null, due to the variation 5<f> being fixed (= 0) at these boundaries, 
as usual. As has been shown in previous works [Izus et ah, 1998; Bouzat & Wio, 1998], it also fulfills the Lyapunov 
characteristic Jj.F[<^] < (it is worth here commenting that, as a matter of fact, this conditions is only valid in a 
weak noise limit). 

Exploiting the so called Hopf-Colc transformation we now define a new field, h(x,t) that, as indicated before, corre- 
sponds to an interface height, 

2v 

h(x, t) = —— ln<f>(x, t), (5) 
A 

with the inverse 

4>{x,t) = e ^( 2 <*). 

As 4>{x,t) > 0, h(x,i) is always well defined. The transformed equation reads 

d t h(x, t) - vV 2 h + ± (Vh) 2 + ^e-%f(h) + £(z, t), (6) 

where f(h) = f(4>) = /(e^ 7 ' 1 )- Now, in order to reduce our general result to the one shown in Eq. ([5]), we assume 

f(<p(x,t)) — acf>(x,t) (with a some constant). In this way Eq. |2]) becomes Eq. (fT]), with a = kt;K and a = (^;) 2 £- 
However, the noise term that in the original equation (Eq. (|2"]l) has a multiplicative character, in the transformed 
equation (Eq. (JTJ) becomes additive. 
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If we now apply the same transformation to the NETLP indicated in Eq.Q, restricting ourselves to f(<j>(x,i)) 
a<p(x,t), we obtain 



T[h]= I e ^.*)-l -K +'-(Vh{x,t)Y dx. (7) 



A 



A 
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It is easy to prove that this functional fulfills both, the relation 



as well as the Lyapunov characteristic -SrJ-[h] < 0, where the function T[h] is given by 



or 

2 



T[h] = 



e 



Hence we have a free energy-like functional from where the KPZ kinetic equation can be obtained through 
functional derivation. Clearly, the contribution to the variation coming from the boundaries is again null. 
It is worth to consider once more Eq. ©, but now assuming f(<f>(x,t)) = a</>(x,i) — &^>(x,i) 3 , i.e. we include a limiting 
(or saturating) term. In such a case, the reaction-diffusion equation corresponds to a form of the so called Schlogl 
model [Mikhailov, 1990; Wio, 1994]. The "potential" contribution to the NETLP in Eq. © has the form 

^(M) 2 - ^OM) 4 ) dx. 

Applying again the indicated Hopf-Cole transformation, in Eq. ([T]) a new associated term arises, having the form 
— jev h ( x,t \ with b = ^7- The new kinetic equation corresponds to a form of the so called bounded-KPZ [Grinstein 
et at, 1996; Tu et at, 1997; de los Santos et al, 2007]. Clearly, we will also have an extra term in the associated 
NETLP (Eq. ([7])). However, here we restrict ourselves to the case 6 = 0, and only analyze the more "usual" form of 
the KPZ equation indicated by Eq. |T]). 



2.2 Some Properties 

Let us now check some of the properties assumed previously for such a functional. According to the analysis of global 
shift invariance in [Hentschel, 1994], it is easy to see that the relations indicated by Eq. (9) in the indicated paper 
arc fulfilled. That is, we can prove that if I is an arbitrary (constant) shift 

F[h + t\ = K[t\F[h] 

T[h + l] = Kil}- 1 ^], (9) 

with K[l] = ei l = (^) 2 T[l}-\ 

To prove other conjectures also indicated in [Hentschel, 1994], we introduce the free energy-like density !F\h, V/i], 
defined through 

F[h] = J dxT[h,Vh]. 

Hence, F[h,Vh] = ±e^ h(s ^[-K + | (Vh(x, t)f]. The relations we refer are 

F[h,Vh] = e s/l J : 1 [(V/i) 2 ] 

T[h,Vh] = e- s/ T![(V/i) 2 ], (10) 

and according to the form of J-[h,Wh], it is clear that the first relation above results obviously true, while for the 
second relation we have that T[h, Vft] = e^'^^Po, where L = 1, and s = ~, as T[h] is not function of Vh. In 
addition, it can be also proved that the indicated NETLP is also invariant under the nonlinear Galilei transformation 
that, as discussed in [Fogedby, 2006], are fulfilled by the KPZ equation. 

It is here adequate to make a warning. We have found that, from the indicated free energy-like functional for the KPZ 
kinetic equation and by a functional derivative, we can obtain a form that resembles a (relaxation) model A according 
to the classification in [Hohenberg & Halperin, 1977]. In one hand, in the standard "model A" it is known that the 
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dynamics can be seen as a superposition of modes that decay exponentially towards a steady state. In addition, it is 
also known that the time-dependent correlations obey certain constraints such as positivity. On the other hand, in the 
KPZ problem it is known that the relaxation of perturbations decay in a stretched exponential manner [Schwartz & 
Edwards, 2002; Edwards & Schwartz, 2002; Colaiori & Moore, 2002; Katzav & Schwartz, 2004]. Hence, even though 
Eq. ([8]) looks similar to "model A" , its behavior is far from trivial, and we have no a-priori intuition to what its 
dynamic could be. Clearly, this is a point to have in mind when suggesting any kind of Antsatz for the temporal 
behavior. 

3. Probability Distribution Function for KPZ 
3.1 General Aspects 

The knowledge of the NETLP indicated in the previous section, allows us to readily write the asymptotic long time 
probability distribution (pdf). We start writing the (functional) Fokker-Planck equation (FPE) associated to Eq. (fT|), 
that reads 

d. f S 

at' 

where, in order to focus in the most relevant aspects and to simplify, we adopted K Q = 0. 

In Chap. 6 of [Barabasi & Stanley, 1995] as well as in Sect. 3.5 of [Halpin-Healy & Zhang, 1995] the form of the 
probability distribution function (pdf) for the one dimensional KPZ equation was discussed. The form of such a pdf 
for an arbitrary dimension, so far, is not known. However, here we show the general form of such a pdf. 
The knowledge of such a NETLP for the KPZ equation allows us to readily write the asymptotic long time probability 
distribution function (pdf), valid for (worth to be remarked) any dimension, which (due to the "diagonal" character 
of T[h]) is given by 

r as [h(x)\ ~ expj-iy dx J* ix) d^rM^ 7 ^ 

r h{x) 



V[h(x,t)} = f dx^{-[vV 2 h(x,t) + ^(Vh(x,t)) 2 ]r{h(x,t)}+e^T[h(x,t)}}-, (11 



6il> J 

exp{-i J dx (r[h]T[h] - dV ^r^M) } 



«W-i^U. (12, 



where, we reiterate, we assumed K = 0, and h re f is an arbitrary reference state. The second line results by using 
functional methods (see for instance [Hanggi, 1985]). The third line shows a nice structure, where we can identify 
a contribution, with a Gaussian dependence on the slope, plus a "correction" term proportional to A. The above 
indicated result is the valid (albeit "formal") solution for arbitrary dimension irrespective of boundary conditions, 
while for the one-dimensional case, and the adequate (periodic) boundary conditions, it is possible to show that this 
pdf reduces to the well known Gaussian result [Barabasi & Stanley, 1995; Halpin-Healy & Zhang, 1995]. 
Clearly, a relevant point is related to the interpretation of the integral over the function ip(x) in Eq. (|12p . In order 
to simplify we consider the 1-d case, and focus only on the first term. Using functional techniques, we have 

J href 



dx 2Ty ' 2J J h ^5ip{x)\J V dx' J I 2 J \dx 



where we have used that |^fey = S(x — x'). Another way to grasp it is via a discrete representation 
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that, as only the I — j and I = j — 1 terms survive, yields 



3 J 



U ^[hj+i-h,] 2 . (14) 
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In both cases we obtain the known result. 

The interpretation for the second (KPZ) term is analogous to the one in the example indicated above. However, 
for this case the situation is a little more delicate. Several recent papers have discussed different alternatives to the 
discrete form of such a contribution in the KPZ equation in 1-d [Newman & Bray, 1996; Lam & Shin, 1998; Ma, Jiang 
& Yang, 2007]. We will not come into these details here, that will be deeply discussed in [Revelli & Wio, 2008]. 

3.2 Nonequilibrium Potential for KPZ 

The last line in Eq. (fT2"]) defines the functional 

r h(x) 



&[h] = - J dx(Vhf - - J dx J d^ (\7^) 2 . (15) 
The variation of this functional give us 



This functional also fulfills the Lyapunov condition J^$>[h] = — ^ sh(s\) ) — ^- 
Hence, such a functional could be identified as the nonequilibrium potential [Graham, 1987; Graham & Tel, 1990; Wio 
et al., 2002] for the KPZ case. More, even though it could be also identified with a Hamiltonian for KPZ, it is worth 
to remark that it differs from the form indicated in Eq. (3.4) of [Halpin-Healy & Zhang, 1995]. 

The last functional form (albeit so far only formal) could, in principle, allow us to exploit several known techniques 
[Graham, 1987; Graham & Tel, 1990; Wio et al., 2002]. It also shows that the claim by previous authors indicated at 
the Introduction is not true. 



3.3 NEP's application: A simple example 

Here we discuss a simple example in order to show the possibilities that offers the knowledge of such a NETLP. For 
this example we analyze a sightly different situation than the one studied in Eq. |T]) and its associated NEP in Eq. 
(fT5|) . We only consider a spatial quenched noise ( "disorder" ) instead of the spatial-temporal noise considered so far. 
The equation associated to such a problem is 

^h{x, t) = vV 2 h{x, t) + ^ (V/i(.x, t)) 2 + K + ${x), (17) 

where, as in previous studies [Nattermann & Renz, 1989; Krug & Halpin-Healy, 1993; Ramasco, Lopez & Rodriguez, 
2006; Szendro, Lopez & Rodriguez, 2007], ft{x) is a quenched, Gaussian distributed, noise. For this case we have that 

with the same form of J-[h] as in Eq. ([7]), but K a replaced by K a + 

The indicated studies have shown that the front profile presents a triangular structure [Ramasco, Lopez & Rodriguez, 
2006; Szendro, Lopez & Rodriguez, 2007] and, clearly, it is of relevance to determine the slope of such structures. 
Using the known form of the NETLP, we can minimize this free-energy-like functional and obtain, in the 1-d case, 

1 /3 

that such a slope is a — f^r) , a value that agrees quite well with the numerical evaluations. 
4. Other Kinetic Equations: Non Locality 



We can go still further and look for the possibility of deriving a NETLP for more general forms of kinetic equations. 
To this end, let us assume that we have the following non-local reaction-diffusion equation [Wio et at, 2002] 



0_ 

at 



<p($,t)= v V 2 (/)(x,t) + a<p(x,t) - f3 / dx'G(x,x')(i)(x',t) + 0(x,t)r)(x\t), 



(19) 



where, as discussed in [Bouzat & Wio, 1999; von Hacften et al., 2004; von Haeften & Wio, 2007], the kernel G(x,x') 
could be of a very general character, and (3 is the interaction intensity. It was shown that the form of the associated 
NETLP is 



?[<!>] = I -2#*> *) + 2 ( V( ^' *)) + £ / dx ' *) x ') ^ *) 



dx. 



(20) 
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As we have done before, using the Hopf-Cole transformation we obtain a generalized-non-local form of the KPZ 
equation 



^h(x, t) = uV 2 h{x, t) + ^ (Vh(x, t)f +K - fi e -^ h{s ^ \ dx'G(x, x~')e^ h(x '•*) + i). 



(21) 



Even though the nonlocal contribution indicated above differs from those discussed in [Mukherji & Bhattacharjee, 
1997; Katzav, 2003], it is clear that such form of a nonlocal term is also of great interest. Repeating the previous 
procedure we find the associated NETLP 



dxe$ h ^ 



A 

2v' 



-K n 



{Vh(x,t)f + p e ~^ h ^ / dx'G(x,x')e^ h ^'^ 



(22) 



At this stage it is required to make some assumptions about the kernel. 

We assume that the nonlocal kernel has translational invariance, that is G(x, x') = G(x-x'). Also, that it is of (very) 
"short" range, that allows us to expand it as 



G(x~x') = J2 A ^ Si2n) ^ 



(23) 



with x — x') — V™,5(x — x') 1 and symmetry properties taken into account. Exploiting this form of the kernel, we 
arrive to the following contributions in Eq. (f2~T|) 



(3 / dx'G(x - x')e* 



A 



2v 

3 



(V/i) 2 V 2 h 



2(g) V 2 (V^-(£) (^h)+^^h 



2v 



2v 



■A* 



(24) 



where the last term indicates contributions of order n > 3 (2n = 6). The parameter (3 could (in principle) be positive 
or negative, indicating an inhibitor or an activator role for the nonlocal interaction term, respectively. 
These contributions have the same form of those ones that arose in several previous works, where scaling properties, 
symmetry arguments, etc, have been used to discuss the possible contributions to a general form of the kinetic equation 
[Hentschel, 1994; Linz et at, 2000; Lopez et ai, 2005]. Clearly, the different contributions that arose in Eq. (j2~4"|) 
are tightly related to several of other previously studied equations, like the Kuramoto-Sivashinsky [Sivashinsky, 1977; 
Kuramoto, 1978], the Sun-Guo-Grant equation [Sun et ai, 1989], and others [Hentschel, 1994]. 



5. Density Dependent Surface Tension 



We can also go further in another direction. Let us consider the case of density dependent diffusion (surface tension). 
Following the same original approach, we start considering the following reaction-diffusion equation with a density 
dependent diffusion for a field 4>(x,t), as studied in [von Hacftcn et ai, 2000; Tessone & Wio, 2007] 

Btftx, t) = V + fW + <f>{x, t)t}(x, t). (25) 

As before, f(4>) is a general nonlinear function, that we collapse to f((f>) = a<f)(x,t). As in Eq. ([2]), r/(x,t) is a delta 
correlated white noise of intensity a. We assume v{4>) — u g((f)), that -using the Hopf-Cole transformation- yields 
u(h) = v g{h). 

The transformed equation reads 

d t h(x, t) = V (v{h)Vh) + ^ (Vhf + K + Z(x, t), (26) 

with X(h) = \ g(h) and a = ^-K a . This equation has a KPZ-like form, with both a density dependent surface tension 
term and a density dependent nonlinear coupling, with both nonlinearities having the same functional dependence. 
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In order to see the associated NETLP, let us remember its form for the original reaction-diffusion equation (Eq. (T25l) ) , 
it reads 



m= J^- jy^')f^')d^ + dx, (27) 



that fulfills [von Haeften et al, 2000] 



d . . I 5T\<j>] , , , 



Applying once more the Hopf-Cole transformation, we obtain 



:F\U\ / ,/r (JjLj |_y due^i>(u)K Q + ^ (v(h)e^Vhyy (28) 
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From this NETLP we can prove that 



where £(x, i) is a white noise of intensity e (with a = (^7~) e )? an d t ne function f [/i] is given by 

-wr- 

As before, if we have f((f>(x,t)) = a<j>(x,i) — b<fi(x,t) 3 , we get a generalized form of the "bounded KPZ" equation. 



6. Conclusions 

The previous results indicates that, for a very general form of a KPZ-like equation, we can device a NETLP a la 
carte. Consider the following equation 

d t h(x, t) = uV 2 h + ^ {\7hf + F{h) + t) (30) 

where F(h) is a general nonlinear function of h. According to the previous results we could readily write the associate 
NETLP that has the form 



F[h] = / { -ir / e^F{u)du + — (e^Vh) } dx. (31) 




Clearly, it fulfills the Lyapunov condition -^!F\h(x, t)] < 0, and also, through functional derivation, allows us to obtain 
the kinetic equation as 

- r ra^ + «*^ (32) 



with T[h] = (^) e-^Hs-t). Hence, it is clear that we can write the NETLP for a very general KPZ-like form, 
independently of the actual form of F(h). 

Summarizing, we have here found the form of the Lyapunov functional or NETLP for the KPZ equation. More, we 
have devised a way to extent the procedure to derive it, and in such a way we were able to derive more general forms, 
including several kinetic equations studied in the literature of interface growing phenomena. Even the case of density 
dependent surface tension, have been discussed. From this NETLP, and through a functional derivative, we have 
obtained either the KPZ as well as other generalized kinetic equations. We have also shown that the NETLP for KPZ 
fulfills global shift properties, as well as other ones anticipated for such an unknown functional. More, we have found 
the exact expression, valid for any dimension, for the probability distribution function, and have commented on the 
result of a simple example that indicate the usefulness of such a functional. 
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As indicated in the literature, dynamic renormalization group techniques, being useful and powerful, in many cases 
only offers incomplete results, having no access to the strong coupling phase [Barabasi & Stanley, 1995; Wiese, 1998]. 
Hence, it is clear the need of alternative ways to analyze the KPZ and related problems, as for instance the self- 
consistent expansion [Katzav & Schwartz, 1999; Katzav, 2003]. The present results open new possibilities of making 
non-perturbational studies for the KPZ problem. For instance, through the analysis of long time mean values of 
h(x,t). In a similar way, it would be possible to obtain correlations, and from them to extract information about 
scaling exponents. Such study will be the subject of forthcoming work. 
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